Non-destructive study of non-equilibrium states of cold, trapped atoms 
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Highly sensitive, non-destructive, real-time spectroscopic determination of the 2D kinetic mo- 
mentum distribution of a cold-atom sample is performed with the three-beam measurement of the 
recoil-induced resonances. The measurements performed with an operating magneto-optical trap 
reveal slow velocity drifts within a stationary atomic cloud and strong anisotropy and asymmetry 
of the non-Maxwellian momentum distribution. The developed method can be easily extended to 
3D. 
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Most experiments with cold, dilute atomic gases em- 
ploy magneto-optical traps (MOT), which yield temper- 
atures in a range of hundreds to a few fj,K. Further traps 
and cooling stages can be applied for reaching the quan- 
tum degeneracy regime. This requires matching of the 
momentum distributions of various traps. Knowledge of 
such distributions is also essential for quantum state di- 
agnostics of the trapped sample. Below, we present reli- 
able method of 2D momentum diagnostics based on the 
so called recoil-induced resonances (RIRs) and apply it 
to the detailed study of non-standard momentum distri- 
butions in a MOT. 

The first unambiguous observation of RIRs was made 
in a ID optical lattice Q] filled with atoms much colder 
than in a standard MOT. RIR signals were also seen in 
optical molasses 0, in a cold atomic beam Q and 
with atoms released from a MOT Q. The influence of 
the recoil effect on the probe absorption and four wave 
mixing spectra has been recently demonstrated in a con- 
tinuously working MOT, i.e. with all light and magnetic 
fields on, in Q. 

In this Letter we present evidence of three different 
kinds of anisotropy of the momentum distribution in an 
operating MOT. The measurements were conducted us- 
ing our three-beam, RIR-based method developed for si- 
multaneous probing of the momentum distribution in two 
perpendicular directions. The method extends the prin- 
ciple of ID thermometry as suggested in Ref. Q. Impor- 
tant feature of our extension is that 2D information is ac- 
quired simultaneously in one measurement. The method 
can be extended to 3D. 

The RIRs result from a stimulated Raman process, 
which couples two kinetic states of free moving atoms 
(Fig. nj. Two laser beams, the pump and the probe, 
of frequencies uj and uj + S, respectively, drive the Ra- 
man transition after which the atoms gain kinetic en- 
ergy Aitkin = fo$ and change momentum p by Ap = 
fiAk = ±2hkei sin6*/2, where k is the modulus of the 
light wave vector, 9 is the angle between the beams, 



and &i is the unit vector perpendicular to the bisector 
of 9. The non-zero amplitude of the considered Raman 
resonance arises from different populations of the given 
kinetic states. When recorded in absorption, the RIR 
shape is proportional to a derivative of the momentum 
distribution dH(pi)/dpi, where pi = p • 0, EJ, 0- 
This direct relation of the RIR signal to n(p) allows con- 
venient and accurate measurement of the kinetic momen- 
tum distribution in a cold atomic sample, provided that 
the distribution is sufficiently narrow. 
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Figure 1: Recoil- induced resonances: (a) geometry of the laser 
beams, (b) atomic kinetic energy and momentum changes due 
to the Raman process coupling two kinetic states. Circles 
symbolize populations II (p) of these states. 

Important advantage of the RIR method is its direc- 
tional selectivity. II(p) is probed in a given direction, 
specified by Ap, i.e. by angle 9 (FigGJi). Hence, apart 
from applications to standard ID velocimetry 0,0], RIRs 
can also be used for studies of a possible momentum dis- 
tribution anisotropy in non-equilibrium states of a cold 
atomic sample. 

Our experiment (Fig. [2j employs a standard vapor- 
loaded MOT @ with 85 Rb atoms. Two extra beams in- 
tersect in the trap center: the pump and the probe. The 
probe beam is directed at a small angle a = 3° to the 
MOT beams (propagating along z), and the pump is at 
9 = 5° to the probe. The probe beam can be detected 
either directly or after retroreflection. The setup with 
retroreflected probe enables the measurement of II(p) si- 
multaneously along two perpendicular directions: _!_, for 
angle 9 between the pump and the nearly co-propagating 
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Figure 2: The layout of the experiment. The pump and probe 
beams intersect in a cloud of cold atoms. With the optional 
mirror we realize a three-beam configuration in which the 
momentum exchange is allowed in two directions: _L and ||. 
Third pair of the MOT beams along x and the MOT coils are 
not shown in this figure. 



probe, and ||, for angle 180° — 9 between the pump and 
the nearly counter-propagating probe. When a and 9 are 
sufficiently small, _L and 1 1 almost coincide with the y and 
z directions, respectively. Both pump and probe beams 
are derived from diode lasers synchronized by injection- 
locking and are blue-detuned from the trapping transi- 
tion 2 S 1/2 (F = 3)- 2 P 3/2 (F' = 4) by A = 2tt-140 MHz 
rs 23.3 T, where T denotes the natural linewidth. Such 
big detuning reduces the perturbation of atoms to a very 
low level (scattering rate oc 1/A 2 ) which is essential for 
non-destructive measurements. Moreover, non-resonant 
pump eliminates overlap of the RIRs and Raman-Zeeman 
resonances Q hence facilitates interpretation of the re- 
sults. Despite large A, the pump and probe beams drive 
the Raman signal with a sufficiently large amplitude and 
signal-to-noise ratio for the pump beam intensities 5- 
35 mW/cm 2 . 

The probe beam is scanned by 5 « ±1 MHz around 
frequency uj of the pump. The probe and pump pho- 
tons induce Raman transitions between atomic kinetic 
states separated by ±h6. Since the polarization of 
the pump and probe beams is chosen to be the same, 
the atoms undergo Raman transitions with Ato_f = 0. 
Hence, the internal atomic state does not change and 
the only states that have to be considered are the ex- 
ternal states associated with the kinetic energy of the 
atomic center-of-mass. The multi-level structure of 85 Rb 
can thus be reduced to a set of independent two-level 
systems which allows straightforward application of the 
basic RIR theory 0, 0, Q|. With the assumption that 
n(p) is the Maxwell-Boltzmann distribution, the RIR 
signal s(8) 0,0] recorded with the retroreflected probe is 
given by two contributions. The narrow one results from 
the Raman process involving the pump and the probe 
beam making small angle 9 and the wide one is for angle 
180° - 9. The signal is 



s(5) oc -A±8exp J - A\\Sexp ^ 2 j 

where, for small 9, gj_ w 2kBk 2 m~ 1 9 2 T± and £ 2 



(S - So)' 



8ksk m~ 7ii. tl and ti are the distribution widths in 
the _L and || directions in the temperature units, A\\ and 
A± are the amplitudes of the corresponding contribu- 
tions, m is the atomic mass, ks is the Boltzmann con- 
stant, and So is the possible frequency shift between the 
_L and || contributions, to be discussed later. 
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Figure 3: Transmission spectrum of a retroreflected probe 
(gray) and the theoretical prediction (dashed) of Eq.(l). In- 
set: the beam setup. The MOT beams have intensity /mot = 
13.8 mW/cm 2 per beam and are detuned from the trapping 
transition by Amot = — 3T, the repumper beam intensity 
/rep = 15 mW/cm 2 , the axial magnetic field gradient d x B = 
12 Gauss/cm, the pump beam intensity / = 33 mW/cm 2 , 
and the probe beam intensity 0.3 mW/cm 2 . 

Typical example of the retroreflected-probe transmis- 
sion spectrum is shown in Fig. It exhibits two distinct 
resonant contributions, predicted by eq. The wide 
contribution is shifted with respect to the narrow one 
by 72.4 kHz, which indicates a 2.8-cm/s average velocity 
component in the || direction. We thus observe an atomic 
drift within a cloud, which as a whole remains station- 
ary. We understand this as a dynamic effect resulting 
from a small difference of the radiation pressures intrin- 
sic to a MOT with retroreflected trapping beams. Indeed, 
the observed shift increases when the imbalance is pur- 
posely increased. Strong imbalance normally produces 
a displacement of the atomic-cloud center of mass. In 
our case this displacement is too weak to be detected by 
standard imaging technique, while the anisotropic atomic 
flow, even one order of magnitude slower than the mean 
thermal velocity, is well measurable with our method. 

As the velocity distributions derived from the signal 
in Fig. El are Gaussian, one can determine the values 
t ± = 172 ± 6 fiK and 7|| = 170 ± 3 fiK. The equality of 
these ts implies thermodynamical equilibrium and allows 
their interpretation as temperature T, despite the slow 
drift. The equilibrium persists for various MOT-light in- 
tensities due to the fact that total intensities of each pair 
of the MOT beams remain the same. The observed nearly 
linear increase of T with the total MOT-beam intensity 
agrees well with previous reports [HI ITlL fl^ . 

The thermodynamics of the system becomes highly 
non-trivial when the trapping light is unevenly dis- 
tributed between the MOT beam pairs. It was predicted 
that for such conditions the width of kinetic momentum 
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distribution shows directional dependence [13j. Using 
the simultaneous measurement of the momentum width 
in two perpendicular directions, we attempted to ob- 
serve such anisotropic non-equilibrium state of the cold- 
atom cloud. For this reason, we changed intensity bal- 
ance between the longitudinal (I z ) and transverse (I x , 
I y ) MOT beam pairs, while keeping the total intensity 
lo = I-x + Iy + Iz constant. We define parameter k as the 



v - (1 - k)7 /2. 
and t± for differ- 



relative intensity of I z , I z = kIq, I x - 
The results of the measurement of T)ig 
ent values of k are depicted in Fig 0] For equal par- 
tition of the trapping intensity (k — 1/3), the widths 
of kinetic momentum distributions are the same, as ex- 
pected. However, when k increases, tii and tj_ follow 
opposite trends, which is evidence of kinetic momentum 
anisotropy in a MOT working MOT and thereby its non- 
equilibrium state. A similar anisotropy was recently ob- 
served also in an optically dense sample 0]. We notice 
that T|| + 2tl, which is the measure of v?< + 2Vj_, is con- 
stant within ±2% over the whole measured range of n. 
The decrease of ri | with the growing k is due to the fact 
that the heating associated with spontaneous emission 
is isotropic, whereas the cooling rate is higher for the 
direction with the increased intensity. The momentum 
anisotropy becomes manifest because the density of the 
atoms is too small to provide efficient thermalization. In- 
deed, simple estimation for typical conditions and Rb-Rb 
collision cross-section CRb— Rb = 3 • 10~ 13 cm 2 j^j yields 
the atomic collision rate below 1 Hz in our trap, while 
the friction coefficient, in frequency units, is in the kHz 
range. 

The theoretical behavior of tm and t± according to 
Refs. [l3l is plotted in Fig. Q] along with the exper- 
imental data. They exhibit similar qualitative depen- 
dence (the decrease of t\\ and the increase of rj_ with 
growing n) , but the existing theory fails to reproduce the 
exact shape of the experimental dependence. This dis- 
crepancy is due to additional mechanism of sub- Dop pler 
cooling, not included in the calculations of Ref . ■ In- 
deed, the increase of I z accompanied by attenuation of l x 
and I y results in efficient quasi-lD cooling scheme in the 
(7+ — optical molasses [JjJ. Evidence of this cooling 
is provided by the values of t\\ falling to 70 jtiK, well be- 
low the Doppler cooling limit of 140 /xK. Importance of 
sub-Doppler cooling for anisotropy of momentum distri- 
bution in cold atomic samples has been previously noted 
in optical molasses [18] . 

In the situation discussed above, the MOT beams were 
carefully aligned which resulted in high stability of the 
trapped-atom cloud, even for the largest departures from 
the equal partition of the trapping light intensity, and al- 
lowed fitting of the RIR signals by eq. (QJ. The sole 
manifestation of the non-equilibrium of the sample was 
the n(pj_) vs. n(p||) anisotropy. The thermodynam- 
ical equilibrium can be altered yet in a different way, 
namely by enhancing imbalance between the counter- 
propagating MOT-beam radiation pressures. Fig. de- 
picts the RIR recorded with a standard ID, two-beam 




Figure 4: Widths of the kinetic momentum distributions 
measured as a function of relative intensity k in the two 
perpendicular directions, t± (hollow circles)and T|| (filled cir- 
cles). For the equilibrium (k = 1/3), /mot = 6.8 mW/cm 2 , 
t± = 100 ± 3 fiK, and t\\ = 103 ± 2 /xK. Other parame- 
ters as in Fig . El Theoretical curves are plotted according to 
Refs. mEl. 
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Figure 5: (a) The probe transmission spectrum in a mis- 
aligned MOT for Jmot = 3.1 mW/cm 2 , Amot = -2.25T 
and other experimental parameters as in Fig EJ As shown 
in inset, the probe-beam is not retroreflected, hence only the 
wide RIR contribution is present, (b) The actual velocity 
distribution II(t)||) (grey thick line) obtained by integrating 
the probe transmission spectrum and the Gaussian reference 
(dotted line), defined in text. 



arrangement applied to the case when the MOT beams 
were slightly misaligned and tuned closer to resonance. 
The ID thermometry was accomplished by replacing the 
optional mirror in Fig. by a photodiode. In this config- 
uration, the pump and the probe make angle 180° — 9 and 
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the recorded signal is proportional only to dH(j>\\)/dp\\. 
Its shape deviates from a derivative of a Gaussian. By 
integrating the signal and scaling to the velocity units, 
the actual distribution of velocity component in the || 
direction, v\\, can be retrieved. Fig. 03 shows such a dis- 
tribution obtained from the experimental signal and the 
idealized Gaussian reference curve of the same area and 
of the width derived from the positions of the minimum 
and maximum of the RIR signal in Fig.|3i. Non-standard 
properties of a stable atomic gas, revealed in our exper- 
iment call for thorough theoretical analysis with proper 
accounting for cooling and heating mechanisms. 

In conclusion, we have developed three-beam spectro- 
scopic method of determining the momentum distribu- 
tions of cold, trapped atoms, based on recoil-induced 
resonances. The method is non-destructive, highly sen- 
sitive and provides multi-dimensional momentum deter- 
mination in a single measurement. Its potential has 
been demonstrated by study of three different momentum 
distributions of atoms in the operating magneto-optical 
trap: (i) thermodynamic equilibrium with well defined 
temperature and Gaussian momentum distribution with 
a slow velocity drift; (ii) non-equilibrium state charac- 
terized by Gaussian distributions with drastically differ- 
ent widths in the longitudinal and transverse directions; 
(iii) non-equilibrium state of non-Gaussian momentum 



distribution along one direction. The result (1 1) q uali- 
tatively confirms theoretical predictions of Ref. and 
indicate need for more refined MOT theory. Our method 
can be particularly useful for studies of anisotropy in op- 
tical molasses [18), 2D MOTs 0, etc. The described 
method can be also straightforwardly applied to 3D case 
by introducing additional pump beam in the xz plane in 
Fig. [2 To avoid overlapping of the RIR signals associated 
with all Ap directions, the frequency of the additional 
pump could be shifted. Such an approach can be used 
for the non-destructive, on-line diagnostics of the atom 
dynamics in a trap carried out independently and simul- 
taneously with other spectroscopic measurements. The 
method should also be applicable to quantum-degenerate 
gases. In fact, the widely used Bragg spectroscopy is 
based on the same principle of momentum and energy 
transfer. Measuring the Bragg-beam transmission in our 
three-beam geometry, rather than imaging BEC can thus 
become a valuable, non-destructive alternative. 
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